On the solubility of Nd 3+ in Y3AI5O12 

D. Klimm a '*, S. Ganschow a , A. Paj^czkowska b , and L. Lipinska 1 

a Institute of Crystal Growth, Max-Born-Str. 2, 12489 Berlin, Germany 
b Institute of Electronic Materials Technology, 01-919 Warsaw, Poland 



Abstract 

Neodymium doped yttrium aluminum garnet powders (Nd:YAG or (Yi-^Ndz^AlsO]^, x = 0.15, 0.25, or 0.30, 
respectively) were prepared by a sol-gel technique. By DTA measurements up to 2000 °C eutectic and liquidus 
temperatures could be determined. Exothermal peaks in the second and subsequent DTA heating runs indicate that 
the crystallized DTA samples are not in equilibrium. The section Y3Al 5 Oi2-Nd3 AI5O12 of the concentration triangle 
Ab03-Nd2 03-Y2 03 is proposed on the basis of thermodynamic calculations that allows to explain the experimental 
results by the balance of metastable phase states in the previously crystallized DTA samples. 
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1. Introduction 

Neodymium doped yttrium aluminum garnet 
(Nd:YAG or Yi_ x Nd x )3Al 5 0i 2 ) is one of the most 
important laser materials. Crystals are available 
from different commercial suppliers with dopant 
concentrations up to 2.5 at% (x < 0.025) [TJ. 
Bakradze et al. [2] (see also reference in [3], Fig. 
04602) studied the subsolidus of the ternary system 
Al203-Nd203-Y203 and found as maximum solu- 
bility for Nd ~ 12 wt-% in the garnet phase, corre- 
sponding to (Yo. 8 2Ndo.i8)3Al 5 Oi2 (x = 0.18). The 
Nd:YAG samples were obtained by cooling from 
PbO/PbF 2 solutions held at 1200 °C. Maximum 
solubility was observed along the "garnet section" 
"G" of the concentration triangle Fig. [TJ 

Recently, some of the present authors reported on 
Nd:YAG powders that were prepared by a sol-gel 
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Fig. 1. Ternary concentration triangle with ternary sections 
G (garnets) and P (perovskites). The solid ternary line limits 
the partial system Al203-Y3AlsOi2-NdA103. 

technique followed by an annealing step at temper- 
atures down to T — 800 °C and found the maxi- 
mum doping concentration (Yo.725Ndo.275)3AlsOi2 
(x = 0.275) [4]. The results were confirmed by X- 
ray and chemical (ICP-AES and SEM-EDX) analy- 
sis. If single crystals were melt-grown from such sol- 
gel powders by the micro-pulling-down technique, 
already x > 0.08 led to the formation of a second 
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(perovskite) phase as the solubility limit of Nd 3+ in 
Y3AI5O12 was exceeded. 

Both reports |2I4| indicate, that the solubility 
limit of Nd in YAG may be large in the subsolidus 
but drops drastically approaching the liquidus 
T > 1900 °C. Unfortunately, the experimental de- 
termination of equilibrium phase relations in the 
subsolidus region is difficult, as the time that is 
needed to reach equilibrium rises exponentially if T 
is lowered. 

The present study will report on DTA measure- 
ments with the same sol-gel samples that were al- 
ready used in the previous study |4| . The results al- 
low to propose the topology of the garnet section 
"G" Y 3 Al50i2-Nd3Al 5 Oi2 of the ternary phase di- 
agram Fig. [TJ 

2. Experimental 

The preparation of nanocrystalline Nd:YAG pow- 
ders (Yi^Nd^AlgOia with x = 0.15,0.25,0.30 
was performed in aqueous solution starting from 
Y2O3 and Nd2C>3 that were dissolved in acetic acid 
and from aluminum nitrate and is described in more 
detail elsewhere [4J. 12 — 18 mg of the apparently 
white Nd:YAG powders were filled in NETZSCH 
standard tungsten DTA crucibles. (The mass of the 
empty crucibles is 500 - 600 mg.) A STA 409C ther- 
mal analyzer with graphite furnace and tungsten 
sample holder (W/Re thermocouples) was used to 
perform the DTA/TG measurements in flowing Ar 
(40 ml/min). The following T program was per- 
formed twice for each sample: 1) From room temper- 
ature to 2000 °C with 15 K/min. 2) Down to 500 °C 
with 20 K/min. 3) Up to 2000 °C with 15 K/min. 4) 
Down to 500 °C with 20 K/min. This means, 4 heat- 
ing/cooling runs were performed in total for each 
sample. It must be noted, that the signal/noise ra- 
tio of DTA signals in the region of very high T < 
2000 °C is large as compared to lower T: 1) A large 
portion of heat is exchanged by radiation and not by 
conduction through the sample holder. 2) The mass 
ratio sample/crucible is worse than for DTA cru- 
cibles that are used for lower T, e.g. from platinum. 

No significant mass loss was found during 
the DTA/TG measurements, indicating that the 
Nd:YAG powders were free of volatiles. After the 
measurements the samples formed solidified melts 
covering the bottom of the W crucibles. The purple 
color of this melt deepened considerably with dop- 
ing level. The DTA curves from the first heating 
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Fig. 2. 1 st (solid), 2 nd (dashed), and 3 rd (dotted) DTA heat- 
ing run of the Nd:YAG samples. (Nd concentration x de- 
noted in the respective graphs.) 

runs of all samples were obviously different from 
the curves of the heating runs 2-4. The cooling 
curves showed strong exothermal peaks at different 
T, indicating undercooling upon crystallization and 
are therefore not well suited for the construction 
of an equilibrium phase diagram. Fig. [2] shows the 
DTA heating curves 1 - 3 of the samples. 

3. Discussion of DTA curves 

A solid that is heated above the fusion point Tf 
melts under consumption of the heat of fusion AHf. 
Accordingly, melting is an endothermal process. If 
the melting phase transformation solid — » liquid pro- 
ceeds through a 2-phase region of the phase diagram, 
the solid melts incongruently. This may be the case 
for mixed crystals (solid solutions, melting starts at 
solidus T SO ] and terminates at liquidus Tn q ) as well 
as for arbitrary compositions within an eutectic sys- 
tem (melting starts at eutectic T eut and terminates 
at Tiiq). All melting events are endothermal. 
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The 15% Nd:YAG sample shows such normal be- 
havior (solid line Fig.[2]top): The first peak with on- 
set T so i = 1764 °C can be attributed to the solidus. 
It is followed by a broad shoulder. The return to 
the baseline indicates termination of the melting at 
the liquidus T Hq = 1898 °C. Contrary, the 25% and 
30% samples do not show any endothermal peaks 
that can be related to the onset of a melting process. 
Instead, broad exothermal effects starting around 
1700 °C can be observed for both samples. 

The 25% Nd:YAG sample shows in the 2 nd and all 
following runs a sharp endothermal effect probably 
related to a eutectic with onset at = 1703 °C. 

For the 30% sample one finds T^ %) = 1700 °C. 
Such eutectic melting is never observed in the 15% 
sample. This behaviour shows, that the composition 
(Yo.85Ndo.i5)3Al50i2 is still within the solubility 
range of Nd3Al 5 0i2 in Y3AI5O12 and the other 
compositions are beyond the solubility limit. Be- 
tween T so i of 15% NdYAG and T eut of the higher 
doped samples some endo- and exothermal effects 
occur until the DTA curve returns to the baseline. 
This return (extrapolated offset) marks the liq- 
uidus that is clearly for the lower doped samples 
T^ %) = 1897°C or T,g 5%) = 1873 °C, respec- 
tively. For (Y . 70 Ndo.3o)3Al 5 12 T^ 0%) = 1860 °C 
was determined by repeated measurements. (Not 
all DTA curves are shown in Fig. El) 

It is surprising that for all compositions the 1 st 
heating curve is greatly different from the following 
ones. Moreover, beginning endothermal peaks do 
sometimes knock over in sharp exothermal effects 
(e.g. 2 nd heating of the 15% sample near 1680 °C). 
Such endo- /exothermal effects during heating show 
bad reproducibility and cannot be explained by 
equilibrium processes. Instead, subsolidus transfor- 
mations must be responsible for exothermal and 
non-reproducible effects, as in the subsolidus range 
the transformation between phases is governed by 
slow diffusion processes and tends thus to pro- 
ceed beyond equilibrium. The following section will 
explain these phenomena on the thermodynamic 
basis. 



the rim system AI2O3-Y2O3. In the opposite rim 
system Al203-Nd203 the perovskite NdAlC>3 is the 
only congruently melting compound. Experimen- 
tal phase diagrams of this system can be found in 
[5f6|7j (see also compilation in [3], Figs. 2342, 6439, 
and 9262), for T > 1300 °C and show in addition 
to NdAlC>3 the /3-alumina type phase AluNdOis 
and the monoclinic (space symmetry group P 2i/c) 
Nd4Al20g that melt both incongruently. 

The "neodymium aluminum garnet" NdaALjO^ 
was, to the authors' knowledge, not yet pre- 
pared. Attempts to produce the phase by the 
same sol-gel technique like the mixed garnets 
(Yi_ x Nda;)3Al50i2, or by crystallization from 
PbO/PbF 2 solutions, failed. Nd 3 Al 5 0i 2 is shown 
as a notional end member in the concentration 
triangle Al 2 03-Nd203-Y 2 03 [2]. Wu and Pelton 
[8] carried out a critical assessment of 15 different 
RE2O3-AI2O3 (RE - rare earth metal) systems, 
among them RE = Nd (but not Y) and calculated 
thermodynamic properties for NdsAlsO]^ and re- 
fined excess enthalpies for the binary melts on 
the basis of a quasi-chemical model. These data 
are incorporated in the data bases coming with 
FactSage [9] and were used for the calculation of 
the Al 2 03-Nd203 rim system in Fig. Q] Following 
this calculation, Nd3Al 5 0i2 should be stable only 
for T < 930 °C and decomposes at higher T to 
NdAlC>3 and a-Al2 03. As the upper stability limit 
of Nd3AlsOi2 is so low one can expect that the 
preparation of this compound will be very difficult, 
if not practically impossible. 

Unfortunately, the very thorough analysis by 
Wu and Pelton [8] did not consider the system 
Y2O3-AI2O3 and well assessed thermodynamic 
data G{T) for Y3AI5O12 are not found in the Fact- 
Sage database. However, the properties of RE-A1 
garnets are known to be very similar, if the ionic 
radii tre are similar too. In octahedral environment 
one finds r Y 3+ = 104.0 pm. The most similar RE 
is holmium with r Ho 3+ = 104.1 pm. The "mixer" 
function of FactSage that is based on a Born-Haber 
cycle allows to calculate data for Y3AI5O12 from 
the notional reaction 



4. Discussion of phase equilibria 

The (Yi-^Ndz^AlsOi^ compositions that were 
investigated in this study can be found on line "G" 
of Fig. [TJ This line starts left at Y 3 Al 5 0i 2 that is 
a congruently melting quasi-binary compound of 



2 Ho 3 Al 5 Oi2 + 3 Y2O3 -> 2 Y 3 Al 5 0i2 + 3 Ho 2 3 

(1) (2) (3) (4) 

where all data except for (3) are known. With (TJ 
one obtains for Y3AI5O12 
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G(T) = -7387342.19 + 3854.78457 T 

-3.7656 x 1CT 3 T 2 + 6326731. 0/T 
-8283.84604 VT - 170451139.0/T 2 
-573.3126 T In T (2) 

c p (T) = 573.3126 + 7.5312 x 1CT 3 T 

-12653462.0/T 2 - 2070.9615/VT 

+ 1.02270684 x 10 9 /!T 3 (3) 

for 298 < T/K < 2213 (G given in J/mol, c p 
in J/(mol K), A#(298) = -7250905.48 J/mol, 
£(298) = 285.461831 J/(mol K)). 

Fig. [3] compares the results of the "mixer" calcula- 
tion with experimental data that were obtained 
by Konings et al. [10] using adiabatic calorimetry 
for T = 5 — 420 K and by drop calorimetry for T = 
470 — 880 K with an experimental error of < 0.35%. 
The theoretical data from |T]) approach the experi- 
mental data around room temperature with a differ- 
ence of 1%. The difference drops with T and reaches 
0.05% at the experimental limit. 
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Fig. 3. Specific heat capacity of Y3AI5O12 calculated from 
JTJ together with experimental data |10j . 

The calculated thermodynamic data G(T) for 
Y3AI5O12 © result at the melting point 1940 °C 
[ITTT2] in a heat of fusion AH ( = 300.9968 k J/mol. 
(Xiao, Derby [13j use for Y3AI5O12 AH ( = 
270.3 kJ/mol without further reference.) Together 
with the data for NdsAlsOi^ from FactSage [9] one 
can calculate the "G" section through the concen- 
tration triangle Fig. [IJ if reasonable estimations are 
done for the 2 mixed phases liquid (I) and garnet 
(g). As excess enthalpies G ex are unknown for both 
phases, ideal mixing and G = Gq + G\d ea \ (Go - 
weighed G for the pure components of the mixed 
phase) was assumed. The influence of neglecting 
G ex was checked for different other RE2O3-AI2O3 



systems where experimental phase diagrams as well 
as Go and G ex are known [319]. It was found, that 
the deviations for liquidus temperatures are typi- 
cally up to T] iq < 150 K; but the topology of the 
whole phase diagram (occurrence of eutectics and 
of 1-, 2- or 3-phase rooms) was never changed. 
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Fig. 4. Calculated section "G" through the ternary phase 
diagram Fig. [T] I - ideal liquid, g - ideal garnet, p - per- 
ovskite NdA103, a - corundum (Q-AI2O3), j3 - /3-alumina 
type AlnNdOig. On the right ordinate the T and phase 
compositions are given, that should be expected from the 
experimental phase diagram Fig. [T] 

The resulting section, that is shown in Fig. [H 
may therefore be regarded as a suitable and self- 
consistent basis for the discussion of the "G" section. 
This section cannot be quasi-binary as the end mem- 
ber Nd3AlsOi2 melts incongruenly and can only be 
understood if the partial system AI2O3-Y3AI5O12- 
NdAlC>3 of the concentration triangle Fig. Q] is con- 
sidered. 

In Fig. |4] the right hand side Nd3Al50i2 corre- 
sponds to the right end of the "G" line in Fig. Q] and 
should basically bear the same information. How- 
ever, Fig. [TJ is based on a thermodynamic assess- 
ment of the whole system Al 2 03-Nd203 by opti- 
mizing G(T) data for all phases, especially G eyi {T) 
for the melt. Unfortunately, the quantity and qual- 
ity of experimental data available for the ternary 
system Al 2 03-Y 3 Al 5 Oi2-NdA103 is not sufficient 
to perform such assessment here too, and ideal mix- 
ing (G ex = 0) was assumed for the melt as well as 
for the garnet phase. Irrespectively of these simpli- 
fications, the correspondence between Figs. [4] and Q] 
is good: For the composition NdsAl50i2 one finds 
the garnet phase stable up to 930 °C as indicated 
by "(g)" near the right ordinate of Fig.[H Upon fur- 
ther heating the garnet decomposes to a mixture of 
the perovskite NdA10 3 and a-Al 2 3 ^ (p+a)'" in 
Fig. [U . At 1393 °C a-Al 2 3 disappears and the (3- 
alumina type AlnNdOig is formed: "(p+P)" in Fig. 
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21 At T = 1713 °C Fig. [TJ proposes disappearence 
of AliiNdOis, resulting in a 2-phase region "(l+p)". 
Fig. S] slightly underestimates the stability of the 
melt with respect to AluNdOis, hence a" (l+p+f3)" 
phase field is found here below " (l+p)". Fortunately, 
this minor deviation contributes only to solid/liquid 
equilibria and no substantial discrepancy in the sub- 
solidus equilibria is evident that are in the focus of 
the present work. 

The (Yi-^Nd^sAlgO^ single phase field " g" 
extends in Fig. [4] from < x < 0.83, but is in 
equilibrium with a liquid phase only for x < 0.12. 
Following this figure, Nd:YAG with doping lev- 
els exceeding « 80% should be stable on the "G" 
section only with additions of a-Al 2 03. For x > 
0.12 (Yi_ x Nda;)3AlgOi2 can only be obtained by 
solid and chemical (sol-gel) reactions. The ternary 
eutectic of the partial system AI2O3-Y3AI5O12- 
NdA10 3 is calculated at r eut , t « 1750 °C. Only 
this ternary eutectic can be observed at constant 
T for different compositions x > 0.28. For 0.12 < 
x < 0.28 the phase boundary "g+p"/"l+g+p" is 
the projection of a binary eutectic groove with 
T eu t,b « 1800 . . . 1750 °C depending on x. In DTA 
heating curves, eutectics are usually marked by 
sharp endothermal peaks. In the present study 
T^ %) = 1703 °C and T^ %) = 1700 °C were mea- 
sured with good reproducibility. The experimental 
values are by « 50 K below the theoretical values, 
but the slightly higher T eut for the 25% sample in- 
dicates, that x — 0.25 is already situated on the 
binary eutectic groove where T eutj b depends on x. 
For large Nd concentrations (depending on T for 
x > 0.12 . . . 0.83) the garnet phase is only stable in 
equilibrium with perovskite, a-A^Oa (corundum), 
or AlnNdOu 08). 

The subsolidus equilibria depicted in Fig. [4] can 
explain the different DTA curves for the first and 
for subsequent heating runs that are shown in Fig. 
[2) If the (Yi_ a ,Nd a; )3Al5 0i2 samples with e.g. x = 
0.30 are heated for the first time, the system starts 
in the u g" phase field and crosses subsequently the 
u g+p", "g+p+P", u l+g+p'\ u l+g'\ and finally the 'T 
fields. The first of the transformations "g" — » "g+p" 
— > u g+p+P" are only transformations between solid 
phases that are diffusion limited and lead to broad 
DTA peaks as the reaction is smeared out over a 
broad T range. If the melt of this composition x = 
0.30 is cooled, the garnet phase crystallizes first with 
enrichment of Nd in the melt, as fc Nd rs 0.15 |14| . 
From Fig. [4]one can estimate fc Nd = a; so i/^ii q ~ 0.12 



for small x. The remaining melt composition ap- 
proaches the eutectic point and crystallizes then on 
the path "l+g+p" -> "g+p+fi" -> "g+p". Indeed, the 
perovskite phase was found in NdYAG fibres that 
were grown in this composition range |4]. However, 
the cooled solid is in a metastable nonequilibrium 
state where Y-rich garnet occurs together with Nd- 
rich garnet and with NdAlC>3 or even AluNdOis- If 
this metastable sample is heated again it tends to 
reach equilibrium by a thermally activated exother- 
mal solid state reaction that is responsible for the 
exothermal DTA peaks. 

5. Summary 

The assumption of the existence of neodymium 
aluminum garnet allows to compile a self-consistent 
YAG-NdAG section of the YsOs-^Os-AkjOs 
phase diagram that explains the measured DTA 
curves for heavily Nd-doped YAG powders. Accord- 
ingly, the solubility of Nd 3+ in YAG ranges up to 
80 at%, i.e. is much higher than supposed so far. 
However, solubility decreases strongly with increas- 
ing temperature. Therefore, and due to distinct 
segregation, the Nd 3+ concentration in crystals 
grown from the melt is limited to < 12 at%. 
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